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We propose and analyze an interface between a topological qubit and a superconducting flux 
qubit. In our scheme, the interaction between Majorana fermions in a topological insulator is 
coherently controlled by a superconducting phase that depends on the quantum state of the flux 
qubit. A controlled phase gate, achieved by pulsing this interaction on and off, can transfer quantum 
information between the topological qubit and the superconducting qubit. 



Introduction Topologically ordered systems are in- 
trinsically robust against local sources of decoherence, 
and therefore hold promise for quantum information pro- 
cessing. There have been many intriguing proposals for 
topological qubits, using spin lattice systems [1 , p+ip 
superconductors [2], and fractional quantum Hall states 
with filling factor 5/2 [3 j. The recently discovered topo- 
logical insulators [4 can also support topologically pro- 
tected qubits [5]. Meanwhile, conventional systems for 
quantum information processing (e.g., ions, spins, pho- 
ton polarizations, superconducting qubits) are steadily 
progressing; recent developments include high fidelity op- 
erations using ions [6] and superconducting qubits [7], 
long-distance entanglement generation using single pho- 
tons [8, 9 , and extremely long coherence times using nu- 
clear spins [TO] . 

Interfaces between topological and conventional quan- 
tum systems have also been considered recently [TTJ [12] . 
Hybrid systems [13] [14] may allow us to combine the 
advantages of conventional qubits (high fidelity readout, 
universal gates, distributed quantum communication and 
computation) with those of topological qubits (robust 
quantum storage, protected gates). In this paper, we 
propose and analyze an interface between a topological 
qubit based on Majorana fermions (MFs) at the surface 
of a topological insulator (TI) [5] and a conventional su- 
perconducting (SC) flux qubit based on a Josephson junc- 
tion device [T5] , The flux qubit has two basis states, for 
which the SC phase of a particular SC island has two 
possible values. In our scheme, this SC phase coherently 
controls the interaction between two MFs on the surface 
of the TI. This coupling between the MFs and the flux 
qubit provides a coherent interface between a topological 
and conventional quantum system, enabling exchange of 
quantum information between the two systems. 

Topological Qubit The topological qubit can be en- 
coded with four Majorana fermion operators 2 3 4 , 

which satisfy the Majorana property jj = ji and 
fermionic anti-commutation relation {7i,7j} = A 
Dirac fermion operator can be constructed from a pair 
of MFs V\- = (ji — ijj ) / y/2 , defining a two dimensional 

Hilbert space labeled by riij = T^T^ = 0, 1. The two 
basis states for the topological qubit, each with an even 
number of Dirac fermions, are |0) t = IO12O34) and 




FIG. 1: (color online). On the surface of TI, patterned SC 
islands can form (a) STIS quantum wire, (b) flux qubit, and 
(c) hybrid system of topological and flux qubits. (a) Two MFs 
(red dots) are localized at two SC tri-j unctions, connected 
by an STIS quantum wire (dashed blue line). The coupling 
between the MFs is controlled by the SC phases <fid = £ and 
4>u = — 7T. (b) A flux qubit consists of four JJs connecting four 
SC islands (a,b,c,d) in series, enclosing an external magnetic 
flux f&o. (c) The hybrid system consists of an STIS wire 
and a flux qubit. The STIS wire (between islands d and u) 
couples the MFs, with coupling strength controlled by the flux 
qubit. The SC phase <j) c can be tuned by a phase-controller 
(not shown), and (j>d = </> c ± 0\ with the choice of ± sign 
depending on the state of the flux qubit. 

|i) topo = II12I34). 

The MFs can be created on the surface of a TI pat- 
terned with s-wave superconductors [5 . Due to the 
proximity effect [TO] , Cooper pairs can tunnel into the 
TI; hence the effective Hamiltonian describing the sur- 
face includes a pairing term, which has the form V = 
Aoe^^j^j + h.c. (where ^j, are electron operators), 
assuming that the chemical potential is close to the Dirac 
point [17]. Here <fi is the SC phase of the island. Each 
MF is localized at an SC vortex that is created by an 
SC tri-junction (i.e., three separated SC islands meet- 
ing at a point, see Fig. [IJa)). The MFs can inter- 
act via a superconductor-TI-superconductor (STIS) wire 
(Fig. [TJa)) that separates the SC islands d and u with 
(fid = e and <p u = — 7r, respectively. For a narrow STIS 
wire with width W <C vp/Ao, the effective Hamiltonian 
is 

H STls = -iv F r x d x + 5 £ r z , (1) 

where vp is the effective fermi velocity, 5 £ = 
Aq cos ((f)d — 4> u ) /2 = —Aq sine/2, and r x,z are Pauli 
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matrices acting on the wire's two zero energy modes [5]. 
As shown in Fig.^a), the STIS wire connects two local- 
ized MFs (indicated by two red dots at the tri-j unctions) 
separated by distance L; these are two of the four MFs 
comprising the topological qubit. The coupling between 
the MFs (denoted as 71 and 72) via the STIS wire can be 
characterized by the Hamiltonian = iE (e) 7172/2, 

with an induced energy splitting E (s) depending on the 
SC phase e. The effective Hamiltonian for the topological 
qubit is 



ttMF 
H 12 



J topo- 



(2) 



where Z topo = (|0> (0| - |1> (l|) topo . 

In Fig. ga), we plot E (e) as a function of a dimen- 
sionless parameter A £ = sin | . For A £ ^> 1 and < 

e < tt/2 0, the energy splitting E (e) w 2 \5 e \ e~ As - 
is negligibly small for localized MFs at the end of the 
wire, as the wavefunctions are proportional to e ~ A £ x / L 
and e ~ A e(L-x)/L _ Q n t ner nanc i ? f or A £ < 1, the 

two MFs are delocalized and E (e) becomes sensitive to 
e. We emphasize that E (s) is a non-linear function of e 
[18], which enables us to switch the coupling on and off. 

Flux Qubit The SC island d can also be part of an SC 
flux qubit (Fig. Qb)), with <fid = s = e° or s 1 depending 
on whether the state of the flux qubit is |0) flux or |l) flux as 
shown in Fig. [2jb,c). Therefore, the Hamiltonian 
couples the flux qubit and the topological qubit. As- 
suming a small phase separation As = e° — e 1 <C 7r/2, 



we can switch off the coupling by tuning e [ 



JQ,1 



to 



satisfy v F /LA <C e ' 1 < tt/2 [5 , so that the MFs 
are localized and uncoupled with negligible energy split- 
ting E (e°) « E (e 1 ) ~ 0. We can also switch on 
the coupling H^ F by adiabatically ramping to the pa- 
rameter regime e 0,1 < vf/LAq to induce a nonnegli- 
gible \E (e°) - E (e 1 )] - A Ae. Because flux qubit 
designs with three Josephson junctions (JJs) [15] [19] 
are not amenable to achieving a small phase separation 
Ae <C 7r/2 (Appendix), we are motivated to modify the 
design of the flux qubit by adding more JJs. 

As shown in Fig. [ljb) , our proposed flux qubit con- 
sists of a loop of four Josephson junctions in series that 
encloses an applied magnetic flux /$ (/ ~ 1/2 and 
^0 = h/2e is the SC flux quantum). The Hamiltonian 
for the flux qubit is 



(3) 



with Josephson potential energy U = 
E;=i 234 (1 — cos 0i), and capacitive charging 

energy T = \ J2i=i 2 3 4^^ 2, ^or tne ^ n is 
the Josephson coupling energy, 0i is the gauge- invariant 
phase difference, C{ is the capacitance, and Vi is the 
voltage across the junction [15] [19]. In addition, there 
are relations satisfied by the phase accumulation around 
the loop ®i + 2/7T = (mod 2tt) and the voltage 
across each junction V{ = (f^) ^ [16] . The parameters 
are chosen as follows: the first two JJs have equal 




FIG. 2: (color online), (a) The energy splitting E (e) (in 
units of AE = vf/L) as a function of A £ = sins/2, 
(b) A contour plot of potential energy U as a function of 
{#1, 02,64} with O3 — 7r — 0\ — O2 — O4. There are two poten- 
tial minima associated with flux qubit states |0) and |1). (c) 
A contour plot of U as a function of {^1,^4} with 0\ — O2 
and O3 = TV — 20i — O4. (d) Marginal probability distribu- 
tions of O4 associated with states |0) (blue solid line) and 
|1) (red dashed line). The parameters are Ej/Ec = 80 and 
{Ej,i/Ej} i=1 ^ 4 = {1, 1, a = 0.8, ft = 10}. 



Josephson coupling energy Ej^ = Ej^ = Ej, the third 
J J has Ej^s = aEj with 0.5 < a < 1, and the fourth 
J J has Ej^ = ftEj with ft ^> 1. For JJs with the same 
thickness but different junction area {^}, Ej^ cx Ai 
and d cx A{. The charging energies can be defined 
as E c ,i E c ,2 = E c = E c ,3 a~ x E c and 

Ec,a = fi~ 1 Ec. For these parameters and / « 1/2, the 
system has two stable states with persistent circulating 
current of opposite sign. We identify the flux qubit basis 
states with the two potential minima |0) flux = 
and J-) flux = |{— 0*}) (modulo 2ir), as illustrated in 
Fig.|gb,c). 

When ft — >• 00, we may neglect the fourth junction 
and this system reduces to the previous flux qubit design 
with three JJs [15J[T9]. For ft ^> 1, there is a small phase 
difference across the fourth J J (Appendix), O4 = ±6\ « 

± where the choice of ± sign depends on the 

direction of the circulating current. We may write O4 = 
Z flux fl|, with Z flux = (|0) (0| - |1) (l|) flux . The fourth JJ 
connects SC islands c and d, and if we fix (j) c relative to <p u 
with a phase-controller [26] , then will be e° = (j) c + 
or e 1 = (j) c — 0% depending on the state of the flux qubit. 



The separation 



As 



1 1 



(4) 



between the two possible values of (j)d becomes small, as 
we desired, when ft is large. 

Aside from this small phase separation, there are also 
quantum fluctuations in O4 due to the finite capacitance. 
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Near its minimum at ± {0*}, the potential energy is ap- 
proximately quadratic; therefore, for f3 ^> 1, the dynam- 
ics of #4 can be well described by a harmonic oscillator 
(HO) Hamiltonian 



H 



HO 



gg4 

2M 4 



J, 4 



(5) 



(at a + 1/2) 



and the canon- 
i (with h = 
cj and 6>4 = 



where the effective mass is M4 = a J- 

ical momentum po 4 satisfies [64, pe^ 
1). We may rewrite i^ HO 
Zflux^J + C (flt + a) / V% where the oscillator frequency 
is uo = ^/8EjEc and the magnitude of quantum fluctua- 

1/4 _ /o n 

Fig. 2[d) shows the proba- 

= e -(#4T#4*) 2 /C 2 



tions is ( 



0" 



-1/2 



bility distribution functions p /i (#4) ~ Cv^ c 
associated with |0)/ and The magnitude of the 

quantum fluctuations £ is comparable to the phase sep- 
aration As; indeed £ ex /3 -1 / 2 may even dominate the 
phase separation Ae cx for large /? (Fig. [3]). [27] 
Therefore, we should consider both the phase separation 
and the quantum fluctuations. 

Hybrid System The Hamiltonian for the hybrid sys- 
tem of topological and flux qubits (Fig. Qc)) is: 



H = H 



HO 



■H™ F = (ata+ 1/2) u-^E(s) Z topo (6) 



where e = </> c +# 4 = </> c + Z flux 6^+C (at + a) />/2. In both 
flux qubit basis states, the oscillator is in its ground state 
with (at a) = 0. To first order in the small parameter 



<C 1, the Hamiltonian becomes 



H = H™-- ((Eo) |0> <0| + (E 1 ) |1) (l|) flux ®Z topo +0 (S 2 

where (E 0/1 ) = J dO^E (<f> c + 4 )po/i (0 4 ). 

Up to a single-qubit rotation, the effective Hamiltonian 
coupling the flux and topological qubits is 



-'flux^topo 



(7) 



with coupling strength g = (Ei) — (Eq) « 
(E (s 1 ) - E ( £ °)) + \ {E" (s 1 ) - E" (e )) C 2 + 
O (£ 3 ) .The first term arises from the phase separation 
and the second term from the quantum fluctuations; 
corrections higher order in £ <C 1 are small. 

Because the energy splitting function E (e) is highly 
non-linear, we may tune <j) c to </> ff such that vf/LAq <C 
= ^ off ± As/ 2 < 7r/2 and switch off the coupling 
g « AoAee-^ off l AoL / 2vF - 0. On the other hand, we 
may adiabatically ramp <fi c to </> on < vf/LAq, which ef- 
fectively switches on the coupling g « AoAe. By adi- 
abatically changing <p c from o fj — >> on — )> o ff with 
j g (t) dt = 7r, we can implement the controlled-phase 
(CPtj) gate between the topological (t) and flux (/) 
qubits. With Hadamard gates Had/, we can achieve 
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FIG. 3: (color online). Comparison between the phase sep- 
aration As oc (dark solid line) and the magnitude of 
quantum fluctuations £ oc /3 -1 ^ 2 (purple dashed line), assum- 
ing Ej/Ec = 80. 



CNOT tJ = Had/ -CP t> / - Had/, which flips the flux qubit 
conditioned on |l) t and can be used for quantum non- 
demolition measurement of the topological qubit [TTJ [20] . 
Furthermore, with Hadamard gates Hadt (implemented 
by exchanging two MFs [H [5]), we can achieve the 
swap operation SWAP t j = (Had; • Had/ • CP t ,/) 3 . Fi- 
nally, with CP^/, Hadt and single-qubit rotations U/, 
we can achieve arbitrary unitary transformations for the 
two-qubit hybrid system of flux and topological qubits 

Imperfections There are four relevant imperfections 
for the coupled system of flux and topological qubits [28] . 
The first imperfection is related to the tunneling between 
and |l)fl ux of the flux qubit, with tunneling rate 

exp ^— yjEjjEc^ . The coupling between flux and 
topological qubits should be strong enough, g ^> t, to 
suppress the undesired tunneling probability r] tU nnei ~ 

(t/gf- 

The next imperfection comes from undesired excita- 
tions of the oscillators. According to the Hamiltonian 
H for the hybrid system, the oscillators may be ex- 
cited via interaction E (0 C + #4) = E (0 C + Zfl ux #|) + 



|0)flux 

t r 



dE > a|+ai 
de S 72 



timated as f] eX cite 

1/4 



The excitation probability can be es- 
Since \m < A , 



\2u de J ' 



C 



(SB*) p 




^SEjEc^ we estimate 



Vexcite 

The third imperfection is due to the finite length of the 
STIS wire, which limits the fidelity for the topological 
qubit itself. When we switch off the coupling between 
the flux and topological qubits by having <j) c = o fj and 
A^ off ^> 1 for the STIS wire, there is an exponentially 
small energy splitting E ~ Aoe _A</, off . 

The last relevant imperfection is associated with the 
excitation modes of the quantum wire, with excitation 
energy E' « vp/L [5 . Occupation of these modes can 
potentially modify the phase separation of the flux qubit. 
Therefore, we need sufficiently low temperature to expo- 
nentially suppress the occupation of these modes by the 
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factor e-^'ABT 

Physical Parameters We may choose the following de- 
sign parameters for the flux qubit: a = 0.8, /3 = 10, 
Ej/E c = 80, and Ej = 200 (2tt) GHz. Both phase sep- 
aration and quantum fluctuations depend sensitively on 
P (see Fig. |3j, with As « 0.16 and C ~ 0.18. Mean- 
while, the flux qubit has plasma oscillation frequency uj « 
60 (27r) GHz, energy barrier AC/ ~ 0.26£?j, tunneling 



rate t ~ l.Sy/EjEc exp 



70 (2tt) 



-0.7 (Ej/E c ) 1/2 
MHz; these parameters only marginally depend on /? 
(Appendix). 

For mesoscopic aluminum junctions with critical cur- 
rent density 500 A/cm 2 , the largest junction (Ej^ — 
j3Ej) has an area of about 1 jim 2 [15] . For the topo- 
logical qubit, it is feasible to achieve the parameters 
A ~ O.lmeV^ 25 (2tt)GHz, v f ~ 10 5 m/s, L ~ 5/ira, 
and T = 20mK. For the interface, the effective cou- 
pling is g ~ AqAe ~ 2(27r)GHz. Therefore, we 
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(assuming <p ^ 



< 



io- d , 

tt/4), 



IS 

have imperfections r] tU nnei r 

e -A0 off ~ e -2O|sin0 off /2| < -j_q 
-3 



and e ~ E ' l kBT < 10~' 6 [29] 

Phase Qubit A similar interface can be constructed 
to couple the SC phase qubit [7] [21] and the topolog- 
ical qubit. A phase qubit is just a J J with a fixed 
DC-current source /. The phase qubit Hamiltonian is 

^phase = T + ^phase^ where j> = and ^phase = 

— 7<l>o0 — Io&o cos 0. The qubit can be encoded in the two 
lowest energy states, |0) phase and |l) phase , with magni- 
tude of quantum fluctuations £q an d Ci ? respectively. The 
coupling strength between phase and topological qubits 
can be estimated as # phase « £"' (e) (d 2 - Co)- 

Notes added It was recently proposed to use 
the Aharonov-Casher (AC) effect for quantum non- 



demolition measurement of a topological qubit [TTJ [12] . 
This proposal, which applies in the parameter regime 
a > 1 where the flux qubit has two possible tunneling 
pathways, exploits the observation that whether two tun- 
neling paths interfere destructively or constructively can 
be controlled by the state of the topological qubit. In 
contrast, our proposal, which applies in the parameter 
regime a < 1 where the flux qubit has only one tunneling 
pathway, exploits the non-linearity of the energy splitting 
E (e) to achieve a controlled-phase coupling between the 
topological and flux qubits. Recently, the related work 
[22] appeared. 



Conclusion We have proposed and analyzed a feasible 
interface between flux and topological qubits. Our pro- 
posal uses a flux qubit design with four JJs, such that the 
two basis states of the qubit have a small phase separa- 
tion Ae on a particular superconducting island, enabling 
us to adiabatically switch on and off the coupling be- 
tween the flux and topological qubits. Such interfaces 
may enable us to store and retrieve quantum informa- 
tion using the topological qubit, to repetitively readout 
the topological qubit with a conventional qubit, or to 
switch between conventional and topological systems for 
various quantum information processing tasks. 
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Appendix A: Flux Qubits 

We describe our design of the SC flux qubit that pre- 
pares an SC island with SC phase e — e° or e 1 depend- 
ing on the flux qubit state |0) flux or |l) flux . The phase 
of the SC island can coherently control the coupling be- 
tween the two Major ana fermions at the end of the STIS 
wire. We would like to have a small phase separation 
Ae = e 1 — e° <C 7r/2, so that we can easily switch off 
the coupling Hfr, F when we do not want to couple the 
Majorana fermions (MFs). 

The previous design of flux qubit with three JJs 
P~5j[IH], however, is not amenable to achieving a small 
phase separation Ae <C 7r/2, because of the following 
reason. The three JJs have Josephson energy Ej^ = 
Ej,2 = Ej and Ej, 3 = aEj [15] [19]. The phase' dif- 
ference across the first junction is = ±cos _1 By 
choosing a = 77 + 1/2 and 77 <C 1, there is a small 
phase separation Ae = 2\0\ ~ 47? 1 / 2 . However, the en- 
ergy barrier for the tunneling is significantly suppressed 
AU = 2a - (2 - Ej « Ar] 2 Ej ~ 4 Ej. The ac- 



tion associated with the tunneling is S « OyJ AU/Ec ~ 
s E j, where the last step uses the prop- 
erty that Ec oc 1/Ej. In order to maintain a similar 
tunneling matrix element between the two potential min- 
ima, it requires the unfavorable scaling Ej oc I/O 3 . For 
example, to achieve = 0.1, we have to increase the 
area of the Josephson junction by ~ 10 3 . Practically, it 
would also be very challenging to have a precise value of 
a = \ + ^- = 0.5025, which fine tunes the phase sepa- 
ration Ae. In contrast, the four-junction design of flux 
qubit has a favorable scaling of Ej^ = [3Ej oc 1/0 and 
there is no fine-tuned parameter. This motivates us to 



redesign the flux qubit with more Josephson junctions. 

As shown in Fig. [4j the flux qubit consists of a loop of 
four JJs in series that encloses an external magnetic flux 
/$o (/ « 1/2 and 3> = h/2e is the SC flux quantum). 
The first two JJs have equal Josephson coupling energy 
Ej y i = Ej y 2 = Ej; the third J J has coupling energy 
^j,3 = &Ej, with 0.5 < a < 1; the coupling in the fourth 
J J is Ej^ = f3Ej, with /3 ^> 1. For JJs with the same 
thickness but different junction area {A,}, Ejj oc Aj 
and Cj ex Aj. The charging energies can be defined as 

e 2 1 

Ec,i = Ec,2 = Ec = Ec,3 = ol~ Eq and Ec^ = 
P~ 1 Ec- Notice that EjjEcj = EjEc is independent of 
j. The system may have two stable states with persistent 
circulating current of opposite sign. 

Here is a summary of the key results: 

1. For j3 — >• 00, we may neglect the fourth JJ 
and reduce the system to the well-studied flux qubit 
with three JJs [15] [19]. For f3 > 1, there is only a 
small phase difference across the fourth junction, with 
O4 = ±01 = Zfl ux #| depending on the sign of the circu- 
lating current (i.e., the state of flux qubit, with Zfl ux = 
(|0) (0| - |1) (l|) flux ). We show that the magnitude of 
the phase difference can be small 0\ w ^taff" 1 ^ • 
As shown in Fig. [4j two SC islands (c and d) are con- 
nected by this junction, if we fix the SC phase C , then 



9d 



-0,1 



± 0% has phase separation Ae = 20% 



aj3 The SC island d can be used to coherently con- 
trol the coupling between the Majorana fermions of the 
STIS wire. 

2. There are quantum fluctuations for the phase of 
the SC island. The magnitude of quantum fluctuations 
depends on {Ecj} and {Ejj}. For (3 ^> 1, the dynamics 
associated with O4 can be characterized by a harmonic 



(a) 



E M — d — 



a f® c 
b 



(b) 
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\ / 

/ \ 






a 




c 








b 







FIG. 4: The design of flux qubit consists of a loop of four 
JJs in series that encloses an external magnetic flux f&o. 
Schematic illustration in terms of (a) JJs and (b) SC islands. 
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oscillator (HO) Hamiltonian 



H 



HO 



Pe 4 
2M A 



Ej,4 



(0 4 



where the effective mass is M4 

ical momentum po 4 satisfies [04,p<9, 
1). We may rewrite i^ HO = (a) a 



1 



and the canon- 
= i (with h = 
1/2) u and (9 4 = 



Zflux^I + C + a ) / V% where the oscillator frequency 
is uj = \/8EjEc and the magnitude of quantum fluctua- 

-1/2 



tions is C = (iff) P 1/2 • We justify that this simple 
model agrees very well with the general model charac- 
terizing the quantum fluctuations for the flux qubit with 
coupled JJs. 

3. Various parameters characterizing the flux qubit 
are also calculated, including the plasma frequen- 
cies {ui} i=12 s, barrier height AC/, and the tunnel- 
ing matrix element t. For example, given parame- 
ters a — 0.8 and f3 = 10, we compute {uoi} w 
(2.8,2.3,1.8)v / S7^ ~ y/EjE^, AU 
" -0.7 (Ej/E c ) 1/2] 



« 0.26-E/j, and 
We notice that 



t « l.Sy/EjEc exp 
these parameters for the flux qubit hardly depend on /? 
when /? > 10, which verifies the intuition that inserting a 
large Josephson junction to the loop has almost no effect 
to the properties of the flux qubit. 

4. We propose two schemes to implement the SC 
phase-controller, which can fix the phase difference be- 
tween an SC island and a big SC reservoir. 

In the following, we provide detailed analysis to justify 
our design of flux qubit with four JJs. First, we give the 
Hamiltonian description for the system. Then, we cal- 
culate the phase separation and quantum fluctuations. 
Next, we numerically obtain various quantities such as 
plasma frequencies, barrier height, and tunneling matrix 
element. Our numerical calculation also verifies our es- 
timates on phase separation and magnitude of quantum 
fluctuations. After that, we propose two implementa- 
tions of the SC phase-controller. Finally, we derive the 
energy splitting function E (s) that is highly non-linear 
in terms of e. 



Appendix B: Hamiltonian for Flux Qubit 

The Hamiltonian for a flux qubit consisting of four JJs 
in series is 



^flux = T + n 



with the Josephson potential 

U = E ^(l-cos^), 

j=l,2,3,4 

and the capacitive charging energy 



T = 



(Bl) 
(B2) 

(B3) 



Here for the j-th Josephson junction, Ejj is the Joseph- 
son coupling, 0j is the gauge- invariant phase difference, 
Cj is the capacitance, and Vj is the voltage across the 
junction. Suppose that all the junctions have the same 



thickness, we have Ejj oc Aj and Ec, : 



2Cj 



oc A, 



where Aj is the area of the junction and Cj is the junc- 
tion capacitance [16 j. The quantity EjjEcj — EjEq 
does not depend on j. 

The Josephson potential is constraint by the phase re- 
lation 



^<9 i + 2/7r = 0(mod27r), 



(B4) 



This phase relation removes one degree of freedom. We 
may introduce a 3- vector = (01,02,03) to character- 
ize the system with four JJs. (Note that we only need 
to choose three independent phases for 3- vector, e.g., 
= (0i,02? #4) is also a valid choice.) For / = 1/2, 
the Josephson potential is 

U = -E Jtl cos 0i - £j, 2 cos 2 - E Jj3 cos 3 (B5) 
+ E JA cos (0i+ 2 + 3 ). 

By appropriately choosing the parameter {Ejj}, there 
are only two minima for the Josephson energy at 
±{0*} = ±{et,6* 2 ,6$,6l} (see Fig. gib). We may 
identify the two levels of the flux qubit as |0) flux = 
|{0*}) and |l) flux = |{-0*}). Thus, we may denoted 
the two potential minima as {Zfl ux 0*}, with Zfl ux = 
(|0)(0|-|l)(l|) flux . 

The capacitive charging energy can be regarded as the 
kinetic energy associated with the dynamics of 0. This 
is because the voltage across the junction is given by the 
Josephson voltage-phase relation Vj = (f^) 0j [16] and 



the time derivatives 



J ~ \2irJ "3 

|0j| obey the constraint 0i 



03 + 04 = (derived from Eq.(B4)). Thus, we may write 
the capacitive charging energy as 



T = 



1 





2 V2tt 



i,j = l,2,3 



1 

= - T • M • 
2 



= -f .w- 1 .^ 



where the capacitive matrix is 



the effective mass tensor is 



M = 



$0 
2tt 



J = l, 2,3,4 



and the canonical momentum is 



p = M • 0. 



(B6) 

(B7) 
(B8) 

(B9) 
(BIO) 

(Bll) 
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FIG. 5: Contour plots of Josephson energy as a function of (a) {^1,^2,^3} (with 64 = tt — 61 — O2 — 63). (b) {^1,^2,^4} (with 
6> 3 = 7T - 0i - 6 2 - 4 ). (c) {0i, 04 } (with 3 = tt - 0i - 2 and 4 = 0) (d,e) {0i, 4 } (with 0i = 2 and 3 = tt - 20i - 4 ). 
(f) Marginal probability distributions of 64 associated with states |0) flux (blue solid line) and |l) flux (red dashed line). The 
parameters are Ej/Ec — 80 and {Ej^/Ej} i=1 2 3 4 = {1, l,a = 0.8,/? = 10}. 



Therefore, we have reduced the problem to the canonical 
model of the quantum system with Hamiltonian 



H = -f rVT 1 - fi- 



ll 



(B12) 



where the operators satisfy the commutation relation 

[0j,Pk] = iSjk- 

Based on this model, we obtain the phase separation 
and the magnitude of quantum fluctuations in the next 
two sections. 



A 



From the phase relation . 6j 
may solve for A. 

For our system with four JJs, we can also calculate 
0| by series expansion with respect to the small param- 
eter For /3 — > 00, we have the zeroth order expan- 
sion A<°) = ^f^, 0< o) = 0i O) = sin" 1 A(°) = cos- 1 
6^ 0) = 7T-2 cos -1 ^ and of } = 0. Then, to the first order 
of P~\ we have 6^ = sin" 1 ^ w 
the phase difference for the fourth junction is 



— . Therefore, 



Appendix C: Phase Separation between the 
Potential Minima 



We now calculate the potential minimum {0*} by 
introducing a Lagr ange variable A associated with the 
phase relation (Eq.(p4j)). We study the function 



- EjjcosQj 

j = l, 2,3,4 



ME' 



(CI) 



The first derivatives all vanish at the extreme point: 



Ejjsme* = A. 



(C2) 



2a/3 



(C3) 



As plotted in Fig. [6b, the numerically obtained quantity 



1 

Ae 



1 

20* 



(X p. 



Appendix D: Models for Quantum Fluctuations of 
the SC Phase 

In this section, we consider two models for quantum 
fluctuations across the JJs. The quantum fluctuations 
across the JJs ar e due to the finite mass matrix for the 
Hamiltonian (Eq.( |B12 )), which is proportional to the ca- 
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FIG. 6: Parameters for the flux qubit as a function of /3. 
(a,b,c) Phase separation Ae (dark solid line) and the magni- 
tude of quantum fluctuations £ (purple dashed line) from nu- 
merical calculation. The theoretical prediction from Eq. ( D3 ) 



(gray solid line in (c)) agrees very well with the numerical 
calculation. (d,e,f) The two ground states of the flux qubit 
can be coupled by intra-cell tunneling (red solid lines) and 
inter-cell tunneling (blue dashed lines) [151 [T9] , characterized 
by the action S, potential barrier AU, and tunneling rate t. 
Panels (a) and (f) assume Ej/Ec — 80. 



pacitance matrix. In the following, we first provide a 
simple, intuitive model to characterize the quantum fluc- 
tuations associated with O4 for j3 ^> 1. Then, we consider 
the general model of multiple coupled harmonic oscilla- 
tors, and obtain the formula for the magnitude of quan- 
tum fluctuations associated with a projected degree of 
freedom. We find very good agreement between the two 
models when f3 ^> 1, which justifies the simple model. 



1. Simple Model — One Harmonic Oscillator 

Given very large /?, we may neglect the higher order 
couplings to other JJs and consider the reduced Hamil- 
tonian for the fourth JJ 

H JA = X -C±Vl + E JA (1 - cos (0 4 - Z flux ^)) , (Dl) 

with displaced potential minimum at Zfl ux 0|. We then 
use the harmonic approximation and obtain 

hHO = it, + {9i ~ Zflu ^) 2 ' (D2) 

1 2 

where M4 = 8Ec - and Ec,4 = 9777. The oscillator has 



frequency uo = \]%Ej^Ec^ = \^8EjEc and characteris- 



tic length 

<=(lff =(t)" v " 2 



(D3) 



Here ( is also the magnitude of quantum fluctuations. 



2. General Model — Coupled Harmonic Oscillators 

We now consider the general model of multiple coupled 
harmonic oscillators, and obtain the formula for the mag- 
nitude of quantum fluctuations associated with a pro- 
jected degree of freedom. 

Close to the potential minimum {Zfl ux #*}, we may 
expand the potential function to the second order of x = 



U = U min + i ^ K ij X i X j + { xS 
i,j=l,2,3 



with 



d 2 U 



13 dOidOj 



(D4) 



(D5) 



The effective Hamiltonian around the minimum describes 
a system of coupled Harmonic oscillators: 



oscillator 



i,j=l,2,3 



-p T • M" 1 • jT+ -x T • K • x 
2 F F 2 



(D6) 
(D7) 

where we have used the definition p = M • = M • ~x . 

To solve for the system of coupled oscillators, we per- 
form the following transformation to make the mass ma- 
trix/tensor isotropic. For real and symmetric mass ma- 
trix (Mij = Mji), there is an orthogonal transformation 
V\ (with V^ 1 = Vi) that diagonalizes the mass matrix 



ViMV{ 



(D8) 



where A^- = \5ij is the diagonal matrix. The eigenvalue 
Xi is effective mass along the i-th principle axis. By in- 
verting both sides, we have ViM _1 V 1 T = A -1 . For di 

agonal matrix, we may also define (A^ 1 / 2 )^ 
Introducing the transformation 



x' = A 1/2 Vi -x 
if = A-V 2 V 1 -p 



(D9) 
(D10) 



we have 



^oscillator = \$ T ' P + \x' T • K • X (Dll) 
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where 



K =A~ 1 ' 2 (V^V?) A" 1 / 2 . 



(D12) 



Finally, we diagonalize the real symmetric matrix K via 
orthogonal transformation V 2 



V 2 KV 2 T = n 



(D13) 



where Q,^ — ojfSij is the diagonal matrix. Overall, the 
position and momentum transform as 

y = V 2 x' = V 2 A 1/2 V 1 -x 



q = V 2 f? = V 2 A~ 1 / 2 V 1 -p 



(D14a) 
(D14b) 



The eigenvalue uJ 2 is the square of the z-th oscillator fre- 
quency (also called plasma frequency). Given parameters 
a = 0.8 and f3 = 10, we calculate the plasma frequencies 
{oJi} i=1 23 = {2.8, 2.3, 1.8} \JEjEc- We may also vary 
the parameter /3, and observe that the plasma frequencies 
only depend very weakly for f3 ^> 1. 

Note that position and momentum have different 
transformations V 2 h 1 / 2 V\ and V^A -1 / 2 "^. Furthermore, 
these transformations are not orthogonal transforma- 
tions. However, as long as the transformations preserve 
the commutation relation [yj,qk] — [ x j,Pk] — ^jfc, we 
can still perform quantization over the transformed co- 
ordinate. 

Following this procedure, we can numerically compute 
the magnitude of quantum fluctuations of O4 as detailed 
below. 



3. Quantum Fluctuations in SC Phase 

We now quantize the phase difference across the fourth 

junction #4. Near the potential minimu m at {#*}, we 

perform the transformation of Eq. (D14) along with the 

It " 

and qi 



quantization iji 



we obtain the 



V2 Hl V2i 

Hamiltonian for three uncoupled harmonic oscillators 



oscillator 



huji ( a\c 
=1,2,3 ^ 



(D15) 



with eigenfrequencies of {wj} i=123 . Each oscillatory 
mode may induce quantum fluctuations in #4, with char- 
acteristic length scale Q. The operator form of 64 can be 
written as 



= Z 



E 

=1,2,3 



aj- 



V2 



(D16) 



We calculate the values of d as the following. In the y- 
coordinate, the characteristic displacement vector for the 

i-th mode is l\ y ^ = \p^ef\ with unit vector along 

the i-th. direction. We may transform this back to the 

x-coordinate, % x) = (V^A 1 / 2 ^)" 1 • Note that 



-f x ) ~ix) 

l\ is no longer orthogonal. From Q ; , we can obtain the 



E 



fc=l,2,3 



characteristic fluctuating scale of Q 

The magnitude of quantum fluctuations of O4 can be com- 
puted as 



1/2 



(D17) 



=1,2,3 



for independent fluctuations from the three uncoupled 
harmonic oscillators. 

As plotted in Fig. |6fc , the numerically obtained value 
for £ (using Eq. (D17)) agrees very w ell w ith the predic- 
tion from the simple model (using Eq. (D3)), wh ich s cales 
as /3 -1 / 2 . Therefore, the simple model of Eq. (D3) pro- 



vides a reliable description to characterize the dynamics 
associated with 64. 



Appendix E: Tunneling Matrix Element 

1. WKB method 

We use the WKB method to estimate the tunneling 
matrix element [19] . The action associated with the path- 
way (r) from 6 (0) = a to (1) = b is 



Je a 



2(U-E)Vd0 T -M-d6, 



(El) 



and the tunneling matrix element can be estimated as 



-S/h 



2tt 



(E2) 



as c/b 

require that the final point, 



The phase space has period 2tt for all three directions. 
We may introduce the unit cell with volume (27r) 3 and 
three basis vectors a\ = 2tt (1,0,0), a 2 = 2tt (0,1,0), 
and as — 2tt (0,0,1). (Note that we may choose the 
shape of the unit cell for our convenience.) Regarding 
the tunneling pathway, we may choose the initial point 
6* = — while the choice for final point is not unique 
^2 i n i a i for integers {n^}. However, we may 
be the one that has the 
minimum action from the initial point, and we choose the 
unit cell so that it includes the minimum action pathway 
that connects #* and 0%. After this procedure, the intra- 
cell tunneling has the minimum action, compared to all 
inter-cell tunneling pathways. 



2. Pathway with minimum action 

We should use the pathway (r) with extreme action 
for the WKB method, i.e., 



SS 
50 (r) 



0. 
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We obtain these extreme pathways using the following 
approach. First, we discretize the integration 



S 



N 



(E3) 



Ej « 200 GHz and Ej/E c ~ 80, which corresponds 
to the first two junctions Ej^ = Ej$ = Ej. Since 
Ejj oc Aj, the fourth junction Ej^ = j3Ej should have 
an area of approximately lxl jam 2 to achieve f3 ~ 10. 



with xo = 6 a and x n = ®b • Then we calculate the deriva- 
tives with respect to fi 



Axi 



f 



Xo 



\Ax 



2 

Xi -\- Xi 



\Ax 



-j 



\Ax i+1 \ v/ 



Xi+l 



(E4) 



For extreme pathway, these derivatives should vanish. 
If we start with an non-extreme pathway with non- 
vanishing derivatives, we can update the pathway so 
that the updated pathway becomes closer to the ex- 
treme pathway. By repeating the update procedure many 
times, we will obtain a pathway that is very close to the 
extreme pathway. Since we know in advance that we are 
looking for the pathway that gives the minimum action, 
we apply the following update rules for the kth update: 



rf +1 >= r -f)_ e 



55 



(E5) 



where e determines the evolution rate and i 
1, 2, • • • , N — 1. For sufficiently small evolution rate 



s[{rf +1) }]-S[{fj k) X] = fSS 



J 



which ensures continuous reduction of the action. 



3. Tunneling rates and f3 Dependence 



O(e 2 )<0, 
(E6) 



Appendix F: Phase-Controllers 

For the STIS quantum wire, we would like to fix 
the phase difference between two disconnected SC is- 
lands. For example, we would like have <\>\ — <p u = 7r/2, 
4>r — = tt/2, <j) c — 4>u — 0c + as shown in Fig. [Ike. 
The idea is to connect the two SC islands via a phase- 
controller. The phase-controller has two large SC islands 
with a controllable phase difference 7. Using large SC is- 
lands for phase-controllers reduces quantum fluctuations 
in the SC phase. 

In this section, we consider two approaches to building 
a phase-controller using Josephson junctions (JJs). The 
phase difference 7 between the two large SC islands can 
be controlled by either the external magnetic flux 7 = 
7 (3>cc) or the electric current 7 = 7 (/), as detailed below. 



(a) 




(b) 



L C=b 




FIG. 7: The phase-controllers can establish desired phase dif- 
ference 7 between SC islands a and b. For example, (a) the 
flux phase-controller with external magnetic flux and (b) 
the current phase-controller with external electric current /. 



This algorithm gives us the correct pathway that lo- 
cally minimize the action between neighboring potential 
wells. We find that the pathways are essentially the 
straight lines connecting the different minima, with no 
significant difference in terms of the action values. For 
example, given parameter s a = 0. 8 and /3 = 10, the intra- 
cell action is Si n w 0.7hy/ 'Ej / Ec and the smallest inter- 
cell action is S out « lAhyjEj/E C - For Ej/E c « 80, 
we will have t 2 /h w exp [(Si - S 2 ) /ft] - 10~ 3 <C 1. The 
barrier height for intra-cell tunneling is AU = 0.25£^j. 

When /3 —> 00, all quantities (potential minimum po- 
sition, plasma frequencies, action for tunneling, energy 
barrier, and tunneling matrix element) reduce to the case 
with three JJs. As illustrated in Fig. [6j the deviation 
scales as j3~ l . For j3 > 10, the perturbation from f3 is 
very small. 

For practical parameters of mesoscopic aluminum junc- 
tions with critical current density 500 A/cm 2 [15] [19], a 
junction with an area of A = 0.2 x 0.4 \im 2 can achieve 



1. Flux phase-controller 

The flux phase- controller uses an rf SQUID loop that 
is interrupted by a single Josephson junction (Fig. [7^1). 
We may change the external magnetic flux enclosed by 
the loop, to induce the desired phase difference between 
the two SC islands, a and b. The rf SQUID loop has 
inductance L, and the J J has Josephson coupling energy 
Ej and capacitive charging energy Ec (= e 2 /2C). The 
Hamiltonian for the flux phase-controller is 

H = U + T, (Fl) 

where the potential energy is 

U =^T (**-P-y) +£j(1-cos 7 ) (F2) 
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with <& x for the external flux enclosed by the loop and 7 
for the gauge invariant phase difference across the junc- 
tion, and the capacitive charging energy is 



1/2 

( 2e/ fie S7 ) • ^he q uan tum fluctuations of 7 around 
the potential minimum has characteristic scale 



T = hjV 2 



(F3) 



$^^r and the effective mass 

Z7T at 



with the voltage V 

m efS = C{^) 2 . 

We find that the potential minimum satisfies the con- 
dition 



$0 
2tt 



7" 



<I>o 



(F4) 



This expression can be used to determine 7 as a function 
of $> x . For L — >• 0, we have 7 = — 27rf^ L . For sufficiently 



small L (satisfying L < 



4ir 2 Ej 



) , this relation is still single 



valued for all 7. Therefore, we can deter ministically con- 
trol the gauge invariant phase difference 7 by applying 
an appropriate <& x . 

We then consider the perturbation around the poten- 
tial minimum and obtain the plasma frequency cu p ~ 

(LC) for L <C 47T 2° Ej • The quantum fluctuations of 
7 around the potential minimum has characteristic scale 



Cflu 



2>A 



1/4 



(F5) 



2. Current phase-controller 

Alternatively, we may also control the phase differences 
via the external current / through a Josephson junction 
(Fig. [7^). The Hamiltonian for such a current phase- 
controller is 



H = U + T, 
where the potential energy is 

Z7 = -/^7 + £j(1-cos 7 ) 

Z7T 



(F6) 



(F7) 



with 7 for the gauge invariant phase difference across the 
junction, and the capacitive charging energy is 



T = -CV 2 
2 



(F8) 



with the voltage V = §^^~ and the effective mass 

As long as / < I c = ^Ej, we have the potential 
minimum at 



7 = sin 1 1/I c . 



(F9) 



We then consider the perturbation around the poten- 
tial minimum and obtain the plasma frequency uj p « 



m eff uj p 



8 E c 
cos 7* Ej 



1/4 



(F10) 



3. Comparison between flux and current 
phase-controllers 

Both flux and current phase-controllers enable us to 
reliably induce phase difference between two SC islands. 
We compare the two phase-controllers in the following 
aspects: 

(1) Tunable phase rage: A single flux phase-controller 
can create all desired phase differences in the range 
[— 7r, 7r]. In contrast, a single current phase-controller can 
only create phase difference in the range (— 7r/2, 7r/2). 
However, by using two or more current phase-controllers 
in series, it is possible to create arbitrary phase differ- 
ences. 

(2) Quantum fluctuations: Both controllers have simi- 
lar scaling for the quantum fluctuations, which scales as 
(E c /E LyJ ) 1/4 . In order to minimize the quantum fluctu- 
ations, we may use the J J loops with a small inductance 
(i.e., El > Ej) for flux phase-controller, while we may 
use the JJ loops with large Josephson coupling energy 
Ej for current phase-controller. Note that for current 
phase-controller the quantum fluctuations becomes un- 
favorably large when 7 w 7r/2, which can be overcome 
by using two or more controllers in series to reduce the 
fluctuations. 



4. Parameters 

We may estimate the quantum fluctuations for prac- 
tical devices. According to the experimental parameters 
of the large Josephson junctions [21]: the charging en- 
ergy E c = e 2 /2C = 0.15 mK and the Josephson cou- 
pling energy Ej = I c ^o/2tt = 500 K, which gives us 
Ccurrent * (8 * ^w^) 1/4 = 0.04. It is also possible 
to build a SQUID loop with very small inductive en- 
ergy E L = $o/2£ = 645 K [25 , and we can obtain 
Cflux « 2y^* (^iiF) 1/4 ~ 0.08. By further increas- 
ing the junction area, we may further decrease E c and 
increase Ej, which should give us more reduced quan- 
tum fluctuations from the phase-controller. 



Appendix G: Brief Derivation for Energy Splitting 

E(e) 

We now briefly derive the energy splitting function 
E(e), which is a highly non-linear function of e. The 
derivation mostly follows Ref. [5]. 



12 



We start with the effective Hamiltonian [5] 

H STls = -iv F r x d x + 5 £ t z , (Gl) 

with S £ = — Aq sin£/2 (differed from [5] by a minus sign, 
due to a slightly different assignment of SC phases). The 
Hamiltonian can be written as 

H s ™ {k) = VpkT * + 8eT * = ( VFk ^ (G2) 

where k is the wave vector in the quantum wire. 
The eigen-energies for H STIS (k) are E^ z (S e ,k) = 
±^S 2 ^v 2 F k 2 . For a finite STIS quantum wire with 
length L, it can only support a discretized set of wave 
vectors satisfying the boundary condition [5 

tan/cL = = kL/A e , (G3) 

with a dimensionless parameter 



For given A e , we may solve Eq. ( |G3| ) and obtain a set 
of solutions kL = f n (A e ) with index n = 0, 1, 2, • • • for 
different bands. The function f n (y) is just the inverse 
function of y = xj tan (x) associated with the nth invert- 
ible domain. 

For the lowest band, we have 

E (e) = E+ (S e , ko) = AE^A* + / 2 (A e ), (G5) 

with AE = vp/L. Note that kL = /o (A e ) is purely 
imaginary for A £ G (1, oo), and it is real for A £ G (— oo, 1). 
Physically, imaginary kL corresponds to localized MFs at 
the ends of the quantum wire, and real kL indicates de- 
localized MFs. Those higher bands (with n > 1) are as- 
sociated with the excitation modes of the quantum wire, 
with excitation energy at least AE [5]. 
In summary, we have 

^ = V^ + ^tAe), (G6) 



A e = -5 e L/v F . 



(G4) 



which is plotted in Fig[2^i. 



